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It has been proposed, that the chiral continuum limit of 2-fiavour QCD with Wilson 
fermions is brought about by a phase in which flavour and parity symmetry are 
broken spontaneously at finite lattice spacing. At finite temperature this phase 
should retract from the weak coupling limit to form 5 cusps. This scenario is 
studied with tree level Symanzik improved actions for both gauge and fermion 
fields on lattices of size 8^ X 4 and 12^ x 24 X 4. 



1 Introduction and Motivation 

The study of the finite temperature phase diagram of 2-flavour QCD with 
Wilson fermions has revealed a rather intricate picture. This picture is Lased 
on the idea of spontaneous breakdown of parity and fla]zpur symmetry El and 
has been investigated analytically as well as numericallyH. The main features 
of this phase diagram are: 

• the critical line Kc{P) defined by a vanishing pion screening mass for 
finite temporal lattice size marks the phase boundary with a phase of 
spontaneously broken parity and flavour symmetry. 

• for large enough Nt five cusps moving towards weak coupling should 
develop, separating the 5 sets of doublers. 

• the finite temperature phase transition line Kt(/3) should not cross the 
critical line, but am past it towards larger values of the hopping param- 
eter, presumablyH turning back toward weak coupling as one crosses the 
niq = line. 

Because the Wilson term breaks all chiral symmetries of the naive fermion 
action, the critical line above cannot consistently be interpreted as thexhiral 
limit of QCD at any finite lattice spacing. It has recently become clearQ that 
two scenarios are possible. Either there exists a second order phase transition 
to an Aoki phase (of width a^) in which parity and flavour symmetry are broken 
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Figure 1: mq(left) and (right) as a function of 1/k on the 12^ X 24 X 4 lattice. 

and along its phase boundary the pion mass vanishes or the system exhibits 
a first order transition along the line of vanishing quark mass and the pion 
does not becomes massless. The analysis also suggests that which scenario is 
realized can change as one varies the action. One therefore has to check, that 
the phase diagram with improved Wilson fermions does exhibit an Aoki phase. 
As the Aoki phase retracts from the weak coupling limit as the temperature 
increases, it naturally explains the absence of any non-analyticities across the 
niq — line in the high temperature phase. In order to separate the high 
temperature side from the low temperature side of the phase diagram, the 
finite temperature transition line should run past the cusp of the Aoki phase 
and continue toward larger k values. The thermal line of the deconfinement 
phase transition therefore crosses the niq = line and should bounce back 
towards weaker coupling due to a swnmetry under the change of sign of the 
mass term in the continuum theoryH. If the thermal line does not touch the 
tip of the cusp, there would be room for the second scenario mentioned above. 

2 The Simulation 

We have simulated 2 flavours of Wilson fermions on lattices of size 8^ x 4 and 
12^ X 24 X 4. We have employed tree level Symanzik improvement which for 
the fermion sector amounts to adding the so called clover term with coefficient 
one and for the gauge fields to adding a 2 x 1 loop. We have first mapped 
out the phase diagram on the smaller lattice and have then corroborated our 
results on the larger lattice for two /3 values. 
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Figure 2: Polyakov loop as a function of k, open symbols for the smaller lattice, filled 
symbols for the larger lattice. 



3 Quark mass and Pion screening mass 

Our results for the quark and pion screening mass are showi in figure (|l|) . The 
current quark mass is defined via the axial Ward identity Q: 
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For /3 = 2.8 we find some curvature for the quark mass as a function of 
though a linear fit produces a reasonable x^. For /3 = 3.1 we have also explored 
a region of hopping parameters where the quark mass becomes negative. There 
we find a rather peculiar behaviour. The quark mass first decreases as one 
lowers K towards Kc and only rises again very close to Kc- This results in 
different slopes of ruq as a function of 1/k for positive and negative quark 
masses, which is in contrast to simulations with unimproved Wilson fcrmions 
that have shown the same behaviour for positive and negative quark massesH. 
We have also accurately measured the pion screening mass. For /3 — 2.8 
the decrease of the pion mass is consistent with a linear behaviour oc 
1/k — 1/kc down to small values of toJ. This also applies for /3 = 3.1 for k 
values that correspond to positive quark masses. For negative quark masses 
the behaviour is quite different and inconsistent with a linear behaviour for 
the points measured. We will however argue below that for /3 = 3.1 the pion 
mass does not go to zero as the quark mass goes to zero, because one crosses 
the finite temperature transition line before the quark mass becomes zero. 
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Figure 3: Chiral order parameter as a function of the quark mass on the 12^ X 24 X 4 lattice. 



4 Polyakov loop 

Figure (Q) shows our results for the Polyakov loop as a function of k including 
data from both lattices. The vertical lines indicate the position of the extrapo- 
lated Kc from the pion screening mass, except for (3 = 3.0 where the pion norm 
was used. This extrapolated Kc decreases with increasing j3. For (3 — 2.8 all 
our data for the pion mass lie in the confined region. We therefore conclude 
that the pion mass vanishes as the quark mass goes to zero, i.e. for (3 — 2.8 we 
hit the Aoki phase as we increase k. For (3 — 3.0 and 3.1 this is no longer so 
clear. The Polyakov loop already shows a high temperature behaviour where 
the extrapolated pion mass would be small. For [3 = 3.0 the situation is less 
prominent and one could still argue that the pion becomes massless, but since 
the Polyakov loop is in the high temperature phase immediately after one 
crosses Kc, the finite temperature line and the line of vanishing quark mass 
come very close together. For /? = 3.1 it becomes evident, that one crosses the 
finite temperature transition line before the line of vanishing quark mass. The 
pion will therefore not become massless as the quark mass vanishes. Because 
the Polyakov loop continues to rise past the ruq — line, one can exclude that 
the finite temperature line bounces back towards weaker coupling immediately. 



5 Chiral order parameter 

Because of the explicit breaking of chiral symmetry by the Wilson action, one 
has to define a properly subtracted chiral order parameter to obtain the correct 
continuum limit. Using axial Ward identities the order parameter is defined 




4 



as follows 

{i>i^)sub = 2mq- Z ■^{TT{x)n{0)) (2) 

X 

Here Z is a renormalisation factor for which we take its tree level value Z = 
(2k)^. The sum over the pion correlation function is just the pion norm. Our 
results are shown in figure (^). For /3 = 2.8 the data extrapolate to a finite 
intercept at niq = 0. For /3 = 3.f the data show more curvature and we expect 
{ip'4')sub to go to zero for vanishing quark mass. 
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